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I. INTRODUCTION 



In this paper we are going to discuss the paradox, concerning the collective quantization in the SU(3) 
Chiral Quark-Soliton Model (xQSM) 0. In the SU(3) xQSM, in analogy with the Skyrme model |@ - 
one constructs the collective wave functions as the eigenstates of the symmetry generators such as 
flavor Ta (or isospin Ta in the SU(2) version of the model)^, and spin operators Sa @- There exists, 
however, the whole set of generalized spin generators including 54. ..7 and iSg which form the SU(3) algebra. 
Precisely speaking 5*8 is a constant equal to —Nc/2^/3. It is however understood, that it is a constant only 
on a subset of the entire Hilbert space, providing a constraint for the physical states. This phenomenon is 
well known already from the SU(3) Skyrme model, and it has been believed that there was no difference 
between the two models as far as the the quantization procedure was concerned. In this paper we will 
show that this equivalence actually breaks down. 

All soliton models start from a classical, statical field configuration U which is subsequently rotated by 
an SU(2) or SU(3) time-dependent matrix A{t). Leaving the details for the next sections, let us quote one 
of the most celebrated results of the SU(3) Skyrme model, namely the relation between the generalized 
spin operators S (called often right generators) and the flavor ones T {left generators): 

% - -D^p{A) Sp, (1) 

where the summation over (3 is assumed. Indeed have proper commutation rules, provided Sp satisfy 
the SU(3) algebra and the coefficients Da/i are Wigner matrices in the adjoint representation of SU(3) 
(or SU(2)): 

AX^A^ = D^p{A) Xp . (2) 

Eq.(|^) has been derived in the Skyrme model and it was believed to be also true in the xQSM. In this 
paper we show that this relation is violated in the SU(3) xQSM: 

= - 1 + DasiA) Ss~Yl ^"/J^^) ' (3) 

\ h / 13=1 



where are dynamical quantities defined in section III, This is precisely the paradox: one constructs 



wave functions as representations of flavor and generalized spin, however, there exists a relation between 
the two sets of generators, which prevents one of them (in our case flavor) from satisfying the proper 
commutation rules. 

It is relatively easy to argue that the troublesome term ^ -* should be discarded, since our calculation 
was not complete as far as all 1/iVc corrections are concerned. Indeed, one can imagine a number of other 
contributions of this order which we have neglected in deriving Eq.(y). There is, however, one catch in 
this reasoning. Namely 5*8 itself is of the order of N^, whereas all other Sp in Eq.(^) are of the order of 1. 
In fact SU(3) quantization also in the Skyrme model mixes different orders in \/Nc expansion. Equation 
(|^) itself contains terms of different order in iVc, nevertheless without the subleading terms, namely 
iSi...7 , the commutation rules for flavor generators would not be satisfled either. For more complicated 
operators, as the charge operator for example, it is not possible to write an analogue of Eq.(|^) and clearly 
decide which terms should be interpreted as non-leading corrections and therefore discarded, and 

which, although subleading, should be nevertheless retained. 

There are two technical reasons for the breakdown of the relation (|l|). The flrst one is that the SU(3) 
soliton models are not fully SU(3) symmetric. Indeed, the SU(3) soliton is usually obtained by the, so 



^Throughout this paper we assume that SU(3) indices, denoted by Greek characters a, /3 . . . run from 1. . .8. 
The SU(2) subset of these indices denoted by capital letters A, B . . . run from 1 to 3, whereas indices denoted by 
a, b . . . run from 4 to 7. 
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called, isospin embedding of the SU(2) soliton. This means that one of the SU(3) generators, in our case 
Ag, has a non-zero expectation value between the soliton states. Let us note in passing, that because 
there is no such operator in the two flavor case, the SU(2) version of the model does not suffer from 
the quantization paradox 0,^. The second reason consists in the fact that xQSM has internal (quark) 
degrees of freedom which, in contrast to the Skyrme model, give rise to the time-nonlocalities 



while calculating expectation values. We will make all these statements clear in section [II, where a toy 
model formulation of the problem will be presented. Earlier, in section ||, we will briefly recapitulate 



semiclassical quantization procedure and give literature where the details can be found. In section [V we 
will reformulate the problem in terms of the path integral. 

Unfortunately, despite many efforts, theoretically satisfactory solution of the quantization paradox has 
not yet been found. Therefore we would like to propose a more modest approach. We will investigate 
a small soliton limit ]lO| of the physical quantities of interest and require that they have limiting values 
corresponding to the nonrelativistic quark model. This will allow us to identify the troublesome 1/iVc 
terms and analyze their structure in terms of the quark degrees of freedom. Detailed analysis of the 
charge operator, axial couplings and magnetic moments will be carried in section 0. For the purpose of 
the comparison with the data we will subsequently discard these troublesome 1/Nc terms, while retaining 
the other subleading contributions. (Symmetry conserving approach) The results will be given in section 



VI , where the model predictions truncated in this way will be compared with experiment and also with 
the SU(2) version of the model, which is free from the quantization ambiguities. 

Although our procedure is in a sense arbitrary, it provides an insight into the origin of the problem. 
Further support comes from the side of phenomenology, since the truncated observables agree better 
(except for the magnetic moments) with the data than the full ones. Moreover, for the quantities which 
can be calculated also in the two flavor case, there is a striking agreeme nt b etween SU(2) and SU(3) 



results. We will summarize our findings and present conclusions in section VII. 



II. A SHORT GUIDE TO THE CHIRAL QUARK SOLITON MODEL 

The generating functional of QCD involves quark and gluon fields. One can imagine the following 
scenario Q first integrate out gluons. The resulting action would then describe the non- linear and non- 
local many quark interactions. The next step would consists in linearizing this complicated action and 
expressing it in terms of local, color singlet composite fields corresponding to the pseudo-scalar mesons 
coupled in a chirally invariant way to the quark fields. And finally we would have to integrate out quarks 
to end up with a pion (or Tr-K-yy) effective Lagrangian. So we go through a chain of effective actions: 

S'qcd[9, A] -> Scs[q] -> -Scff [g, tt, K, rj] Scs[tt, K, 77]. (4) 

It should be kept in mind that the arrows in Eq. (^ do not indicate a rigorous derivation of one action 
from another but rather educated guesses based mainly on symmetry principles and some physical input. 
As 5cff [(?, TT, if, 77] we will choose a semibosonized action of the non-linear Nambu-Jona-Lasinio (NJL) 
model, which is follows from QCD in the instanton liquid model of the QCD vacuum [|l2|. It is based on 
the following quark- meson Lagrangian (Chiral Quark-Soliton Model): 

L ^ ^{ij^'d^, - m - MU^''')il^ , (5) 

where is defined by = ^^^^ U + ^-^^U^ and M is the constituent quark mass of the order 
350-450 MeV. In what follows we will neglect current quark masses rh. 

Next we will assume that baryons can be described as solitons [|l^ - of the effective action cor- 
responding to the Lagrangian (g). That means that the Goldstone fields described by matrix U are, in 



^See recent review by D.I. Diakonov ||l] 
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a sense, large and fulfill classical equations of motion. To this end we choose a hedgehog Ansatz for U, 
which in the case of 2 flavors takes the following form (see reviews: |9|,p^): 

U2{r) = ex.p {ir ■ fi F{r)) . (6) 

Here F{r) is a profile function, and n — r/r is a unit vector. In SU(3) case the soliton field reads: 



t/3(r) 



U2 
1 



(7) 



Obviously U3 is not completely symmetric in the SU(3) flavor space; namely it commutes with As- 
Next one allows for the rotation in the flavor space: t/(r) — + A{t)U{'r)A^ (t) with angular velocity 



~iTr 



The Lagrangian of the system takes then the following form: 



L = ^P^A (^idt + m • V - /3M C/^' - ^Xa^o^ = ^j^'A D[U, Q] A^ . 



(8) 



(9) 



The detailed path integral derivation of the collective Hamiltonian describing baryon spectrum will be 
given in section where the subtleties connected with the non-local nature of the effective action will 
be discussed. By integrating out the quark fields one arrives at the collective Lagrangian £, which is a 
power series in VL: 



2\/3 



I2 ^a^6 



(10) 



where matrices /i^2 are the inertia tensors in the subspaces 1 ... 3 and 4 ... 7 respectively. Normal proce- 
dure consists in constructing effective Hamiltonian 7i acting in the representation space spanned by wave 
functions, which transform as irreducible representations of two symmetry groups: 1) left multiplication 
of A by an SU(3) matrix and 2) right multiplication of A by an SU(2) matrix. The left symmetry corre- 
sponds to flavor, whereas the right one to spin. In fact one can promote the right symmetry group to the 
full SU(3) and then consider only the subset of the full Hilbert space, for which the right hypercharge 
is equal to —1. This constraint follows from the fact that Jig appears only linearly in the collective 
Lagrangian (^o|). Therefore the baryon wave function of flavor^ B = {YTTz) and spin S = {—1883) can 
be explicitly written in terms of the SU(3) Wigner D^^' functions: 



in) 

BS 



D 



{YTT3)ilS-S3 



(11) 



Note that the lowest SU(3) representations TZ allowed by the constraint Yr = — iVc/3 correspond to octet 
and decuplet for Nc — 3. 

In section IV we will show in more detail how to calculate tensors of inertia and expectation values of 
various operators in the path integral formalism. Before doing that, let us consider a simple quantum- 
mechanical model, which reveals all necessary features needed to state the quantization paradox. 



^Obviously in the SU(2) case index a is confined to 1. . .3 and Gell-Mann matrices A should be replaced by Pauli 
matrices r. 

''Here three numbers in brackets correspond to the SU(3) states of the left (B) or right (S) symmetry group. 
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III. QUANTUM-MECHANICAL TOY MODEL 

Consider^ a relativistic Hamiltonian H[U], corresponding to the Lagrangian (|^), satisfying a single- 
particle eigen-equation: 

H[U]\z)=e,\t). (12) 

In Eq.(|l^) the color index of the quark fields has been suppressed. As usually we assume that the ground 
state consists of all Dirac see levels and one positive energy level, called valence, filled. 
Rotation induces a perturbation into the Hamiltonian (p^): 

H[U] A (^H[U] - ^Xo^Qo^ = AH'[U]A^. (13) 

For further calculations it is convenient to work in the rotating frame: 

|^)=At|^), (14) 

where the eigen-equation reads: 

H[U]-^Xan^y^)^e^\^). (15) 

We will solve Eq. (^5|) perturbatively, noting that the unperturbed spectrum satisfies Eq. ([l^ ) . In the first 
order we get: 

Now in order to calculate the matrix element of any intrinsic operator consisting of a string of 7 
matrices (and possibly some space dependent function) F and a flavor matrix Ac,: 

Ol = KT (17) 

we have first to transform it into the rotating frame: 

Ol^dl = AOlA^ = D^pXp F . (18) 

Operator Oa has to be subsequently sandwiched between the occupied states ( p^ ) and summed over i. In 
this way we get expression for the collective operator Oa which acts in the Hilbert space of the collective 
baryon states given in Eq.(^l]): 

Ol= HDapXpTlm) 

m^occ. 

_ 1 ^ {m\DapXpT\n){n\Xi3ni3\m) + {m\Xf3fli3\n){n\DapXpT\m) ^^^^ 



2 

nGnon — occ 



n 

mGocc. 



Throughout this paper we shall always assume that n runs over non-occupied states, whereas m runs over 
occupied states (including the valence level). Of course the above formula for 0~ requires regularization 
which does not concern us for the moment. 



^We thank W. Broniowski for pointing to us this simple way of stating the problem. 
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Note that in Eq.(p^ matrix Dap{A) and angular velocity appear in two different orders. We are 
not going to commute them because in the process of collective quantization f2^ will be promoted to the 
angular momentum operators, which do not commute with Dap{A). Although it may seem that there is 
some arbitrariness in the way we define 0~ as far as the ordering of Dap{A) and fi^ is concerned, we will 
stick to the natural order dictated by the perturbative expansion of the rotated levels. We shall argue in 
the next section that in the path integral formulation the order of the operators is actually unambiguous. 

Now using decomposition into commutator and anticommutator: 

Dap^P = -[Dap,^f3\ + -{Dap,^(3} , 

npD^p^ -]^[D^p,np] + ]^{D^p,ni3} (20) 



we can rewrite Eq.(]l9|) as follows: 



— {D^p, Up) 2^ -— 



^ m.n ^" ^™ 

where we have explicitly pulled out an overall factor corresponding to the sum over colors. 

Equation ( pi]) is our final result for any observable. Let us first consider collective isospin generators 
in the SU(2) version of the model (with obvious replacement of Gell-Mann matrices A by Pauli matrices 
r). Of course we expect to reproduce Eq.(|^) in this case. For the isospin F = 1/2 and the first term in 
Eq.(|2^) is zero since the SU(2) hedgehog is fully symmetric. For the same reason the last term vanishes 
and we are left with the simple formula: 

{DabM (22) 

where the SU(2) moment of inertia is defined through the following relation: 

iVc ^ {m\Tc\n){n\TB\m) + {m\TB\n){n\Tc\m) ^ ^(+) ^^^^ 



In order to reproduce Eq.(|i|) we have to define spin generators. We shall do this first for SU(3) and 
then restrict the general formula to the SU(2) case. To this end let us observe that the total angular 
momentum Ja is defined ab initio in the rotating frame. Therefore there is no need to transform spin 
generators^ to the rotating frame, as it was done in Eq.([T^). Because of the hedgehog symmetry the 
matrix elements of J a are equal to the matrix elements of —1/2 Xa, which allows us to define generalized 



spin generators corresponding to all eight A matrices. In fact we can simply use Eq.(21) with F — 1/2, 
replacing Dap — > —Sap, which gives: 

c i\ I \^^<=n (m|Aa|n)(n|A^|TO) + (wlA^ |n)(n|Aa|m) 

ba ^ -— 2_^{m\Xa\m) + —ilp ; ■ (24) 



One should remember that Ja, which correspond to the total angular momentum of the soliton, are subsequently 
quantized as spins of the collective states. That is why we interchangeably call them spins or total angular 
momentum operators, which, however, should not cause any confusion. 
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Note that the antisymmetric part (the third term in Eq.(piD) does not contribute here. 

Let us turn to the SU(2) case. Again the first term in Eq.([24|) vanishes and in the second term we 
immediately recognize the moment of inertia defined in Eg. (p3|) . FinaUy, due to the commutation 
rule 



[Sp,Dap] ~ ifpp-fDa-f or [Sb,Dar\ — i^BRcDAC 



(25) 



the order of the operators in Dab and Sb (summed over B) in the anticommutator in Eq.(|22|) does not 
matter and one arrives, as expected, at Eq.(|^). 

What goes wrong in the SU(3) case? Let us first observe that Eq.(p^ separates into two independent 
parts: one corresponding to the subspace of the SU(2) soliton (see Eq.(|^)) and the plane wave sector 
corresponding to the "strange" quark ^ Plane waves with negative energy which arc occupied will be 
denoted by mo, whereas the non-occupied positive energy levels will be denoted as no- 



soliton states 




pi 



ane waves 



T3i 
O 

o 

h 
o 



mo 



T3i 

^ > 

O CD 

O - — I 



FIG. 1. Spectrum of the Dirac Hamiltonian (12). 



Let us first calculate Sa ■ To this end let us observe that the first term in Eq. (|2J) does not vanish since 
As has a non-vanishing expectation value: 



A"^^^(m|A„ 



\m) = — ^ das, 



2^3 



(26) 



provided one subtracts contribution of the vacuum. Next, the sum multiplying fJ^ splits into two inde- 
pendent parts: 



iVc {m\K\n){n\\(3\m) + (m|A^|rt)(n|Aa|TO) (+) (+) 

-J- = OaB l-i + Oab l2 



(27) 



This equation defines l'^^ , which can be reexpressed in terms of the overlaps between the soliton states 
and the plane wave states: 



'^This nomenclature is not entirely correct, since we work in the rotating frame, where all quark flavors are 
mixed. 
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^ (m|no)(no|m) ^ ^ (TOo|n)(n|mo) 1 ^^g) 

These transitions ale illustrated in Fig.l. 

Throughout this paper we use the subscript " 1" in the definitions of the inertia parameters, if only the 
soliton states are involved, while if subscript "2" is used, it is understood that transitions between the 
plane waves and the soliton states also contribute. Superscripts " (±)" are used to denote symmetrized or 
antisymmetrized contributions. Clearly l[ ' = 0, because the soliton is totally symmetric in the SU(2) 
subspace. 

Finally we get the expression for Sa '■ 

Sa - -^^"8 + + ■ (29) 



Now we are ready to calculate flavor generators. Substituting F = 1/2 in Eq.(|is|) we observe that: 

Nc o„iY^ {m\Xp\n){n\Xp\rn)_-jjTi\Xi3\n){n\Xp\rn) 



T n t A\ Q rn n i {m\Xp\n){n\Xp\m) - {m\Xfi\n){n\Xp\m) 

ht = -Uap{A) bj3 —[Dap, lip] — ■ (oU) 




This is precisely the last term in Eq.(|3C|) which causes the problem. Were fi's the classical quantities, 
the commutator would vanish. 

It is easy to convince oneself that the last term in Eq. ( ^0|) is unequal to zero only for the indices p — r 
and P = b , with r, 5 = 4 ... 7: 

Nc {m\Xr\n){n\Xb\m) - {m\Xb\n){n\Xr\m) .2 (_) 

m.n * 

Again, I2 can be rewritten in terms of the overlaps between the plane waves and the soliton states: 

^ (TO|no)(no|m) _ ^ (TOo|n)(n|mo) | ^^^^ 

It is now clear that for the fully symmetric Ansatz ' would vanish. Indeed, in that case the plane 
wave spectrum would have to be replaced by the soliton spectrum, and the two terms in Eq.(|3^) would 
cancel. This is precisely the reason why there is no antisymmetric piece in the SU(2) case {i.e. for p = R 
and p = B , with R,B = 1...3). 

The final expression for the fiavor generators is obtained by contracting / symbols: 

/sbc fbcp — 3 (Jsp, (33) 

yielding: 

r(-) 

T„ = -Dafi {A) Sfi + V37^„8 {A) ^ . (34) 

^1 

A few remarks concerning the last term in Eq.(|3^) are in order. First, this term spoils commutation 
rules for the flavor generators. Moreover it is non-diagonal in the SU(3) representation space. Finally, it 
would vanish for the symmetric spectrum, as clearly seen from Eq.(p2|) and Fig.l. 

Using relation (Q) one can calculate the nucleon electric charge by sandwiching the charge operator 
between the hadronic states of Eq. (0) : 



e^=(7^3 + M+Ur3 + ?')^^, (35) 



2 / 5 V 2/ /2(+) 

where T3 stands for the isospin of the nucleon. 

Finally, let us note that generally for F 7^ 1 there will be antisymmetric contributions in Eq.(^l|) both 
within the soliton states (i.e. for p = i?, /3 = B, with R,B = 1 ... 3) and for the transitions between the 
soliton states and the plane waves. 
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IV. PATH INTEGRAL FORMULATION OF THE PROBLEM 



A. Collective Lagrangian and adiabatic approximation 

In this section we will rederive Eq.(|o|) and Eq.(pTl) within the path integral formalism. As wc shall 
see, the path integral will uniquely determine the order of the Dap functions and angular velocities fi^g. 
To this end let us consider the 3rd action from the chain of Eq. , with the Lagrangian density given by 
Eq.(|^). The pertinent partition function reads: 

j V^'^V^PVU e'l'^"''^ . (36) 

Let us first outline steps needed to derive the collective Lagrangian (^0|). For this purpose one considers 
the nucleon correlation function: 

n^(f,T/2;y,-T/2) = (o|jA,(f,T/2)J^(y,-r/2)|o) , (37) 

where the JN{x,t) and jl^{x,t) are nucleon annihilation and creation operators, respectively. It is 
implicitly assumed that the nucleon consists of Nc rather than of 3 quarks. Although we work in the 
Minkowski space-time we will take the limit iT ^ cxd in order to select the ground state contribution to 
the correlator (|37|): 

lini Un{x, T/2; y, ~ exp{-iTE,i) . (38) 



Within the path integral formahsm Eq.(|37|) reads: 

n^(f,T/2;y,-T/2) ^^J V^^VibVU Jjv(.t,T/2) e'/'^'"'^ jj,{y,-T/2) . 

The integration over the quark fields can be carried out exactly. In order to integrate over the chiral 
meson field U we first introduce a classical solution, the soliton, which is given by a dominant trajectory 
in the large Nc limit. The classical equation of motion for U is derived from the variation of the classical 
energy E^i in Eq.(^8|). In this way the soliton profile function F and the soliton mass M are calculated. 

Next, we take only zero-frequency modes of U around the classical solution, i.e. the rotational modes^. 
Then, the path integral VU is replaced by V^, where are Euler angles corresponding to the rotation 
matrix A{t) = A{£,{t)): 

d'xd^y UN{x,T/2;y,-T/2) ^ j *jv(e(r/2)) e^^SpiogD[c/,«] ^],^(^(_t/2)) , (39) 

where Sp stands for the functional and matrix trace of the logarithm of the Dirac operator D[U, S] defined 
in Eq.(^. The ^'Af(0 ^ nucleon collective wave function, which is a function of the Euler angles ^q. It 
is normalized by: 

ON'N ■ 

It is already clear from Eq.(|39|) that the effective action (which, by the way, corresponds to the last 
one in Eq.(y)) is highly nonlocal. This fact is further visible if one expands SplogD in powers of the 
angular velocity fl up to O(il^). This is pictorially illustrated in Fig. 2 where the quark loop with the full 
propagator is expanded into the powers of fl represented by crosses. 



We neglect here translational mode which would restore translational invariance. 
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• • • 



FIG. 2. Full quark loop (thick line) expanded in powers of angular velocity SI (denoted by crosses). 
Because of the time dependence of fl, higher order terms give a nonlocality in time: 

J d^xtfy nw(x,T/2;y,-T/2) ^ j *]v(C(^/2)) e~*^=' t+ j dt,dt,c^^(t,M) ^^(^(_t/2)) 

The collective nonlocal Lagrangian £NL(ti,^2) should be understood as: 
j dtidt2 Cni.{tiM) j dtidt2{-A"nM Siti-t2) 



(40) 



+ 7;^aitl) I"^{h ~ t2) np{t2) e{tl - t2) 



(41) 



where 



i(e,i-em)t 



TO|AQ|n)(n|A/3|r7i) , 



where states n and m have been already introduced in the previous section and A" is given in Eg. (|2q) . 

In order to recover the local Lagrangian ( [lo| ) we have to make an adiabatic approximation and assume 
that Oq. very weakly depends on time. This amounts to: ^ait2) ~ ^a{ti). HencePi 



J dtidt2 I'^ih - 12) ^^a(^2) oih - 12) ~ 

J dh dt2i''^{h-t2)e{h-t2)^nM 

The integration of the I°'^{ti — ^2) over t2 gives: 

— —t 

j.,^ linr 7 dt/"/^(t) = ^y ^"^'^"'"^^"'^^'"^^ ■ 

»T-,oo J 2 f-; £„-£,„ 



(42) 



(43) 



Note that this integration is finite only if Sm < £n- This condition and the Pauli principle make that only 
transitions from occupied m to nonoccupied n states are possible. 

Needless to say that formula ( ^3| ) contains both symmetric and antisymmetric contributions. Symmetric 
contribution has been defined in Eqs.(p3|, [2^ ), whereas the antisymmetric one is given by (31). However, 
due to the symmetricity of the product ri^rj^ the antisymmetric part does not contribute to the collective 
Lagrangian. 

The nonlocality in time encountered here is relatively harmless. By employing adiabatic approximation 
we recover the local collective Lagrangian (|T^ and the collective quantization can be carried through in 
the normal way. Let us also mention that in many cases one employs a different procedure to calculate 
moments of inertia using the identity for the real part Re logD = 1/2 logD^D. This method, although 
useful 6.17. in proper time regularization, obscures the problem of time nonlocality. In general, however, 
the time nonlocality will bring up some new terms which will be of importance for us here. We will see 
this when we consider nucleon matrix elements of local quark operators. 



Local and nonlocal inertia parameters can be distinguished by the explicit time dependence in the latter case 
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B. Nucleon matrix elements 



Let us consider a local quark operator O defined in analogy with Eg. ([191): 



(44) 



We will calculate now the nucleon matrix element of the collective operator in the path integral 
formalism: 



{N,p'\Oi{ti)\N,p) =^ fd^xd'ye+'P'-^e-'Pyfvi;^ 



(45) 



Jn{S,+T/2) Olih) e^/^'^'^j;(y,-r/2) , 

where the Z' is a normalization factor which gives: 

{N,p'\N,p) = {2TTfS'{p'~p) . 

After integrating over the quark fields and the collective coordinates (rotational zero mode and also 
translational zero mode), we obtain 



dt 



Dc^iiih) ip{x, h - 12) n^{t2) e(ti - 12) 



(46) 



+ n,{t2) 4'' Hx,t, ~ t2) D^pih) e{t2 - h)] } e»/^*^(*) ^n{-T/2) 



Note that time dependence of D matrices, fi's and wave functions 5* at appears through time dependence 
of the Euler angles = (,{t). The densities introduced in Eq.(^6|) are defined by: 



A(^(f) = Nc^{m\x)rXf3{x\m) 



i(en-em)t 



m\x)TXi3{x\n) (n\X~f\m) 



IP\x,t)=^^Y. 



ra\Wn) {n\x)V\i3{x\m) 



(47) 




FIG. 3. Matrix element of the external current denoted by a curly line expanded in powers of angular velocity 
Q, (denoted by crosses). 

Equation ( ^6|) defines clearly the order of the operators fJ-y and Dap- This is pictorially illustrated in 
Fig.||. Only now, following the trick of adiabatic approximation employed already in Eqs.(^ 43), and 
using ( po| ) we arrive at the formula for the matrix element of the operator which has been already 
anticipated in Eq.(pl|). Let us remark that the formulae obtained with the help of Eqs.(|47|) are slightly 
more general - they allow for calculation of momentum dependent matrix elements, which we will need 
to calculate form factors, for example. Integrated over dx with p' = p they reduce to (OT 
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Let us now apply our results to the operators of physical interest. In section III we have already 
calculated the charge operator. Next, let us consider axial vector current. To this end one considers an 
operator: 

O'^Zk = i^^AlKl^i^ ■ (48) 
Finally, we will consider magnetic moment operator: 

^?QK =V^AQ(rX7)^V- (49) 
Here Q corresponds to the electric charge: 

Q + ^As . (50) 

Note that operator (^) depends on space variable r and it involves nontrivial integration over (Px. 

Let us analyze various contributions entering Eq.(^l|) and define the pertinent inertia parameters. From 
general symmetry considerations one gets: 

-i[Do.p,np] (^-^eKB.BX[-^ + fKrt X^-^^ . (51) 

Equation ( pl| ) defines inertia parameters x[^^ which depend explicitly on the Lorentz structure F. They 
can be read off directly from Eq.(|2l|). After some algebra we get: 

OaK = \ Xo - -^JJ- + '^f+y j ^"''^ ^ DanSx 2dpqKDapSq . (52) 

This form of O^j^ can be also used for the singlet operators, i.e. for a = provided one makes a 
replacement: 

-D08 — * V3 otherwise Dop . (53) 
Finally, let us write explicit forms of the inertia parameters both for the axial current [ p^ : 

1^ (m|(g X f)3|7i)(7i|r3|m) , 



m,n 



^(±) ^ ^ I ^ {m\a-f\no){n„\m) ^ ^ {mo\a ■ T\n) {n\mo) I ^^^^ 



and for magnetic moments 
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m[~^ = (m|((f X (t) X f)3\n){n\T3\m) , 



4 £"7^ Ef-fi 

771, n 



^(±) = J ^ (m|75(fx a) ■ T\no){no\m) ^ ^ (mo|75(fx a) ■ T\n){n\mo) I ^^^^ 

Our analysis in this section has been based on the unregularized effective action (|36|). For practical 
calculations one has to employ technically more involved methods in order to find the explicit form of 
the regulators. The relevant formulae can be found in the literature |17 1^. However, not all inertia 



parameters require regularization. In the case of the axial current and the magnetic moment operator 
A'{^2 s-nd ^12^ do not have to be regularized. They follow from the imaginary part of the action in the 
Euclidean formulation of the model and they have counterparts in the Skyrme model jl^ where they are 
related to the Wess-Zumino term. Quantities with superscript "(— )" do not have counterparts in the 
Skyrme model|^. They are related to the non-locality of the effective action. 



V. SMALL SOLITON LIMIT 



We shall discuss now the collective operators for the axial couplings and magnetic moments. Collective, 
in the sense that they have to be sandwiched between the baryon wave functions (pi]). Equation ( |5^ ) 
together with the expression for the charge (given already in Eqs. (|34|j35| )) suffer from the quantization 
dilemma discussed in the introduction. In this section we will study the behavior of various terms in 
Eqs.(|5^) and ( ^5|) in the limit of the small soliton. 

Chiral Quark-Soliton Model interpolates between the naive nonrelativistic quark model and the Skyrme 
model. To this end one studies the model with fixed parameters (constituent quark mass M and cutoff 
parameters) by varying the soliton profile function F(r). If the soliton size tq (defined through F = 
F{r/ro)) decreases, the valence level joins the upper continuum and the sea contribution vanishes. For 
large solitons the valence level sinks into the Dirac sea and the explicit contribution of the valence level 
disappears. In Ref. it was shown that in the SU(2) version of the model the well known quark model 

result gf ^ = (N^ + 2)/3 can be recovered only if the rotational correction proportional to a[ V^i^"* 
taken into account. This is a convincing argument in favor of the rotational corrections which originate 
from the time nonlocality of the quark loop. Moreover, in the limit of large soliton size this correction 
is dying out much faster than the leading contribution. This result is also welcome, since in this limit 
xQSM results should reduce to the ones of the Skyrme model. Indeed, in the Skyrme model this kind of 
rotational corrections does not exist, because the model is based on a local Lagrangian density. 

Here we are going to reverse the logic used in Ref. [|l^ . We will investigate the small soliton size of the 
SU(3) formulae ( p2[ ) in order to identify terms which spoil the agreement with the quark model results. 
We will also demand that both versions of the model SU(2) and SU(3) have the same limit for small 
solitons. These two requirements suffice to identify unambiguously the troublesome terms. As we shall 
see these terms correspond to the quantities Aj of Eq.(|5l|). 

The nonrelativistic quark model (NQM) limit is most naturally defined in terms of a somewhat unphys- 
ical profile function F = 7r6{ro — r) with tq 0. In this limit the energy of the valence level i?vai M 
and the sea contributions to all quantities vanish. Simultaneously we take the nonrelativistic limit, i.e. 



Precisely speaking in the simplest version of the Skyrme model without vector mesons. 
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(±) 



we systematically neglect small components of the spinors | i ). Although the inertia parameters 

are infinite in this limit because of the energy denominators, their ratios to the moments of inertia /{^ 
remain finite. In particular: 



(+) 



'■2 



(56) 



and the nucleon charge (35) has wrong NQM limit as tq — > 0. 

Let us next consider the axial coupling . To this end we have to calculate matrix elements of the 
collective operators involving Wigner D functions in the state of a proton with spin up. The results for 
SU(2) and SU(3) take the following form: 



[SU(2)] = -:lMo- 



l 



A 



r(+) 



5f [SU(3)] = -- 



A 



(-) 



A. 



(-)' 



1 A 



(+) 



30 /(+) 



Next we consider NQM limit of the quantities entering Eq.(p 



A, 







-3 , 



A. 



(+) 



1 , 



A 



(-) 



r(+) 



2 , 



A 



(±) 



7 A. 



(+) 



30 /(+) 



(57) 



(58) 



With these values we obtain 

,(3) 



ffrisu(2)] 



and 



fff [SU(3)] 



where 7/15 for gj^'' [SU(3)] is the limiting value of A2 ' /I. 

Finally we will consider magnetic moments of proton and neutron: 



7_ 
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AiN[SU(2)] 

m[su(3)] 



1 M- 



(+) 



6 



M- 



(-) 



r(+) 



2n 



1 

60 

1 

60 



Mo 



Ml 



(-) 



M. 



(+) 



M. 



M- 



(+) 



(+) 



r(+) 



(+) 



7 M. 



M- 







(+) 



+ 



M. 



(+) 



M 



(+) 



r(+) 



r(+) 



2r. 



(59) 



(60) 



Inertia parameters are given in terms of the space integral K which can be found in the Appendix. In 
the limit of the small soliton size we get: 



M, 







-2K 



M- 



(+) 



(+) 



-K 



M- 



(-) 



(+) 



-K 



M, 



(±) 



r(+) 



-K 



With these values we obtain 

1 



1 



Atjv[SU(2)] ^ -(l±5)if and /iAr[SU(3)] -> -(1 ± 5)i^ + — (1 ± 7)i^ 
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where + refers to the proton and 



45 



(61) 



to the neutron, respectively. The last term in Eq.(|61|) corresponds 
1-2 1^-2 - It can be readily seen from Eq.(|6l|), that the quark model ratio 
fJ-p/fJ-n — —3/2 is recovered only if this term is neglected. Indeed: 



to the hmiting value of M^^ V^^'' 
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where 1/13 corresponds to M2 1 12 ■ 

In this section we have calculated small soliton limit of the collective operators corresponding to charge, 
axial decay constants and magnetic moments. In the SU(2) case the limiting values coincide with the ones 

of the nonrelativistic quark model. In the case of SU(3) they get extra contributions coming from X2 
defined in Eq.(|l]). If we neglect these terms the nonrelativistic quark model values are recovered. Let us 
stress that the agreement with the quark model results cannot be obtained without the contributions of 

X[^\ which originate from the time non-locality of the fermion loop as X2 ^ ■ There reason why X2 "* do 
not vanish, is the isospin embedding of the SU(2) hedgehog in the SU(3) U field given by Eq.(|^). Because 

of this asymmetric embedding the Dirac spectrum ( p^ is asymmetric, and, as a result, X2 ' are non-zero. 
In the case of the charge operator it is relatively easy to argue, as it was done in the introduction, that 
I2 corresponds to the nonleading 1/iVc correction to the coefficient in front of Sg,. However, it is 

not possible to make such a clear distinction for other quantities. Nevertheless, in the following we will 
neglect X2 ^ terms. In that way the agreement with the nonrelativistic quark model will be restored and, 
as we will see in the next section, we will be able to describe the experimental data with much better 
degree of accuracy. 

Let us finally note that even if we discard terms proportional to X^ ^ the SU(3) expressions, in 
comparison to SU(2), contain extra terms, as seen explicitly from Eqs.(|5^) and (|60|). Nevertheless for 
ro ^ the limiting values for both SU(2) and SU(3) do coincide. This is due to the fact that the matrix 
elements of the D functions are much smaller in the case of SU(3) and the extra terms make up for this 
difference, so that finally the SU(2) and SU(3) limiting values are the same. This will be also almost 
exactly true for the realistic self-consistent soliton profile. 

VI. NUMERICAL RESULTS 

The method we use for the numerical computation is based on the sums over quark single-particle levels 
in the background pion field. The system is in a spherical 3-dimensional box where the eigenfunctions of 
the Dirac Hamiltonian in the background pion field can be obtained by numerical diagonalization [ pO| . 
For the ultraviolet regularization of the Dirac sea contribution we employ the proper-time scheme. The 
ultraviolet cutoff A is determined by fitting the pion decay constant and it reads: A ~ 600MeV. The 
soliton profile is found by self-consistent minimization of the static classical energy. All results are for 
m^r = mK = 140 MeV. 




0.0 ^ ' , ' , ' , ' , ' 1 

0.0 0.2 0.4 0.6 0.8 1.0 

QlGeVl 



15 



FIG. 4. Proton electric form factor in the SU(2) (dashed line) and sym. cons. SU(3) (solid line) version of the 
model. The experimental data are from ref. |2l|. 
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FIG. 5. Neutron electric form factor in the SU(2) (dashed line) and sym. cons. SU (3) (solid line) version of 
the model. The experimental data are from ref. |22| denoted by solid circles and ref. p3] denoted by an open 
triangle. 

Before considering SU(3) results within the prescription which we have developed in this work it is 
instructive to review the SU(2) results given in Ref. Fig. ^ shows the proton electric form factor 
calculated with M — 420 MeV. The dashed curve corresponds to the SU(2) result and it follows the 
experimental data with very good accuracy. In Fig. || we plot neutron electric form factor. The dashed 
curve corresponding to the SU(2) version of the model overestimates the data. We will come back to this 
point later. In Fig. ^ we plot the axial form factor of the nucleon. It is also in a good agreement with 
the experiments. In this case the next-to-leading order correction of 1/Nc expansion of the rotational 

zero-mode of soliton, namely A\~'^ in Eq.(|5^), is crucial to reproduce the experimental data |@,||. 
Moreover this correction ensures that in the limit of the small soliton our result agrees with the estima- 
tion of the naive non-relativistic quark model: 5^'' = {Nc + 2)/3 [0. In Figs. |^ and || the magnetic 
form factors of proton and neutron are shown. For the magnetic properties we also take into account the 
next-to- leading order correction of 1/Nc expansion, namely m[ ^ in Eq.(|60|). The theoretical curves 
are multiplied with a factor such that the values at = agree with the corresponding experimental 
magnetic moments. Apparently the g-dependence of the form factors is well reproduced. The theoretical 
magnetic moments deviate noticeably from the experimental ones, as one can see at Table. This unfor- 
tunate feature is in accordance with most, if not all chiral models which notoriously have difficulties to 
obtain proper absolute values of the magnetic moments. 





SU(2) 


sym. cons. SU(3) 


cxp. 


^p[n.m.] 


1.98 


1.81 


2.79 


Mn [n.m.] 


-1.36 


-1.20 


-1.91 
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FIG. 6. Nucleon axial form factor in the SUf 2) (dashed hne) and sym. cons. SU(3) (sohd hne) version of the 
model. The experimental data are from ref. |25[|26|] . 

The solid curves in Figs. |^|| correspond to the SU(3) results truncated according to the prescription 
advocated in section 0. The differences between the results in SU(2) and SU(3) are very small, except 
the case of the neutron electric form factor, Fig. H. 
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FIG. 7. Proton magnetic form factor in the SU(2) (dashed line) and sym. cons. SU(3) (solid line) version of 
the model. The experimental data are from ref. |2l|] . 
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FIG. 8. Neutron magnetic form factor in the SU(2) (dashed line) and sym. cons. SU(3) (solid line) version 
of the model. The experimental data are from ref. denoted by solid circles and ref. |2j] denoted by open 
triangles. 

In the xQSM all quantities consists from two kinds of contributions. One is the valence level con- 
tribution, and the other one corresponds to the Dirac sea levels. The proton electric form factor is 
dominated by the valence level contribution, while the neutron one is dominated by the Dirac sea, which 
is interpreted as contribution of a mesonic cloud. The remarkable difference between SU(2)- and SU(3) 
calculation exists only for the neutron electric form factor and there are two possible explanations: 1) It 
might be due to a systematic error of the model. Because the neutron electric form factor is quite tiny 
compared with the one of the proton, it is quite probable that it is sensitive to the systematic error of 
the model. 2) It might be due to a kaonic excitations of the vacuum. In the SU(2) the vacuum excitation 
(Dirac sea polarization) in the neutron occurs only in the tt" channel. On the other hand, in the SU(3), 
it has another contribution corresponding to the if"*" channel. The excitation suppresses the tt" 
excitation. Because of that, the SU(3) result is below the result in SU(2). One should also stress the 
sensitivity of the magnetic form factors to the tail of the soliton field. 

In this paper we have treated strange degrees of freedom in a purely perturbative way. One possible 
method which assumes strong SU(3) breaking is given by the bound-state approach proposed by Callan 
and Klebanov in the Skyrme model |27| . We plan to investigate the neutron electric form factor with a 
method based on the bound-state approach to the kaon. 

VII. SUMMARY 

The aim of this paper was to discuss in full extent the quantization paradox, which arises, when one 
applies canonical quantization rules to the rotating soliton in the SU(3) Chiral Quark-Soliton Model. 

We have argued, that canonical quantization leads to the two type of terms, called x\ 2 j which appear 
due to the non-commutativity of the SU(3) Wigner functions Dap and the angular velocities fi^, and 
the non-locality of the fermion loop. These terms constitute 1 /Nc corrections to the physical observables 
in question, and, at first sight, it seems that in the consistent treatment either both types should be 
simultaneously retained or neglected. However, as we have shown in full detail, quantities x[ \ which 
involve only the transitions between the states from the SU(2) soliton subspace are theoretically harmless 

and phenomenologically desperately needed, whereas terms X2 ^ , which involve the transitions between 
the soliton states and plane waves lead to the violation of the canonical quantization rules for the flavor 
generators. So far no satisfactory theoretical solution of this paradox has been found. We have proposed 
a semi-phenomenological way, based on the requirement that the xQSM model results should agree with 
the non-relativistic quark model in the limit of the artificially small soliton, to circumvent the apparent 
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contradiction and discard the ^2 terms. The only justification of this procedure, which we can give at 
the moment, is based on the fact, that the SU(3) canonical quantization mixes different orders of 1/Nc 
expansion and it is a priori not clear which terms are only formally subleading, and which constitute 
genuine corrections to the physical quantities. We have discussed this at length in the Introduction 
on the example of the nucleon charge, where such a distinction is easy to make. Unfortunately other 
observables like axial couplings and magnetic moments, which we have discussed in this paper, are much 
more complicated, and the same kind of argument cannot be straightforwardly applied in that case. 
Once one accepts the philosophy that x[~'^ terms should be kept as they come out from the calculations 

and X2 ■* terms should be discarded, a consistent picture of the SU(3) soliton emerges. All collective 
operators obey proper commutation rules, and all observables have proper (i.e. non-relativistic quark 
model) limit for very small solitons. The SU(2) and SU(3) results coincide and the agreement with the 
experiment is, generally speaking, satisfactory. 

We felt obliged to discuss this problem without being able to give sound theoretical solution to it, since 
it plagues SU(3) model calculations published in the literature. We feel that non-expert readers should 
be aware of this problem, with the hope that some of them may find an elegant explanation in a not 
distant future. 
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